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y^Q '. Abstract 

. Deser et al. proposed a combination of the Einstein and Landau-Lifshitz 

C\| ! pseudotensors such that the second derivatives in vacuum are proportional to 

O ' the Bel-Robinson tensor. Stimulated by their work, the present paper dis- 

cuss the gravitational energy-momentum expression which has the same de- 
sirable Bel-Robinson tensor property. We find modifications of the Einstein 
and Landau-Lifshitz pseudotensors that both give the same coefficient of the 
Bel- Robinson tensor in vacuum in holonomic frames. 



0> 



1 Introduction 

o 

The classical pseudotensor is not a tensorial object, it is frame dependent. A de- 
sirable requirement is that to the lowest non-vanishing order in vacuum the energy 
density for a gravitating system should be proportional to the Bel-Robinson tensor 
O . B a/3nv This will assume a covariant and positive energy on all reference frames, 

i.e. the associated energy-momentum vector is future pointing and non- spacelike. We 
find one expression associated with a good, but more complicated boundary condition 
which gives the Bel- Robinson tensor only [2]. However, it is not easy to obtain the 
Bel-Robinson tensor alone and in general it is accompanied by certain other tensors, 
S a pfivi Yapfiv and T a p^ y [3]. 
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According to the standard textbook [I], the Einstein pseudotensor does not give 
a positive gravitational energy in vacuum to second order . Recently Deser et al. 
[3] found the analogous expansion for the Landau-Lifshitz pseudotensor by using the 
similar method, however this too does not have the expected result. Yet these two 
classical pseudotensors are good, giving physical sensible results, inside matter (mass 
density) and at spatial infinity (ADM mass). There comes a natural question can 
they be modified, so that they have the extra nice result of the positive gravitational 
energy in vacuum? The answer is yes. 

Moreover, the present paper gives the coefficient of the positive gravitational en- 
ergy for some suitable pseudotensors expressions in holonomic frames, namely the 
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combination of the Einstein and Landau-Lifshitz pseudotensors [3], one of the Chen- 
Nester four quasilocal expressions [2] and Papapetrou pseudotensor [5]. 



2 Ingredients 

2.1 Riemann normal coordinate 

In general the metric tensor g a p, using the Taylor series expansion, can be expanded 
as 

g«p(x) = g aP {0) + {d^g aP ){0)x» + ~(<9^ a/3 )(0)x^ + . . . (1) 
At the origin in Riemann normal coordinates 

g a p(0) = Vap, d^g a/3 (0) = 0, (2) 

Hence ^ 

5W?0) = Vaf3 ~ -R a ^K^x K + 0(x 3 ). (4) 

2.2 The tensors B, S, Y, T and their significance 

In order to extract the physical meaning of the tensors B al 3^ u , T a p^ V) S a p^ v and Y a p^ V) 
one can use the analog of the "electric" E^ and "magnetic" parts in vacuum of 
the Weyl tensor, 

Eij = Coioj, Hij = *Coioj, (5) 
where is the Weyl conformal tensor and *C a p^ u is its dual, 

* C a /3fiu = -V—9 € af3\aC Xa (6) 

In a simple form using the Riemann tensor in vacuum 

E a b = RoaOb, H a b = -^R^amn^b^ ■ (7) 

Certain commonly occurring quadratic combinations of the Riemann tensor compo- 
nents in terms of the electric E a b and magnetic H a b parts in vacuum are 

D r> a b 771 rpab p r> abc n tt rjab r> -pabcd at? rpab fo\ 

-rtOaOfe-Ko — &ab& > J^Oabc^O — ^n a btl , ti a bcdti — ^ab^ ■ [o) 
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In particular, the Riemann squared tensor can be written in terms of the electric and 
magnetic parts as 

p TDOtfiiw a p pOaOb i a p pOaftc , p ryabcd 

— ^-KOaOfe-K + 4H 0abc K + ^abcd^ 

= 8(E ab E ab ~H ab H ab ), (9) 

where Greek letter means 0,1,2,3 and Latin stands for 1,2,3. The Bel-Robinson 
tensor B a p^ u and tensor S a p^ v are defined as follows 

'■— RaXiujRp v* + RaXuaR/3 // ~ -QafiQ 'fiuRprXaR^ °\ (10) 

o 

S a f3^iv '■= RatiXa Rf3u Xa + RavXoRfi^ + -^jafiQiivRprXaR^^ ■ (H) 

For our future analysis, define the tensors Y a ^ v and T a/3fll/ [3] 

:= RaXfiaRi^ J 7 + RaXfia Rjy X if , (12) 
: = — — 9al39^RxapTR XapT - (13) 

The physical observable of the energy- momentum for the tensors B a fs^ u , T a/3fiu , 
and Yafi^v are the spatial traces 

B^oi 1 = (E ab E ab + H ab H ab , 2e c ab E ad H d b ), (14) 

T^i = (E ab E ab -H ab H ab ,0), (15) 

S^oi = -lOT^o/, (16) 

Y^oi = 2(E ab E ab ,e c ab E ad H d b ). (17) 

Note that B^qi has a form similar to the Maxwell energy-momentum density and in 
particular -Boo/ > 0. There is one relation which needs our attention throughout the 
present text 

This means that it is not necessary to obtain the Bel- Robinson tensor B^ 1 for the 
positive energy requirement, a combination of the tensors Y^oi 1 with T^oi or equiva- 
lently Y^i with S^qi can suffice. 



Buqi — S U Qi + Y u0 i + 9T u oi . (li 



3 The interior, ADM and gravitational energy 

Note three physical quantities of interest if one considers a massive object in general 
relativity. They are the interior mass density, the ADM mass [S] at the spatial 
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infinity and the gravitational field energy-momentum in vacuum. In order to study the 
gravitational energy, Einstein proposed the classical pseudotensor t a ^ which follows 
from the superpotential Uj^. Unfortunately the superpotential is not uniquely 
defined, for example 

tot = d v V^ v \ (19) 
but one can introduce a new pseudotensor such as 

t a » = t a » + d u Ua [M , (20) 

which is likewise conserved may seem that there is no special way to study the grav- 
itational energy. However, the interior and ADM mass provided some guidelines or 
restriction, so that one can have some sort of physical energy-momentum components. 
Recalling the Freud superpotential, one may consider the generalization 

u a w = ^ {fci(<S£rV - *«rV) + h(5^\ - 5^\) + k 3 (r^ a - r^)} , 

(21) 

where k±, k 2 and k 3 are the extra added constants. Inside matter at the origin in 
Riemann normal coordinates 



^{(h + 2k 2 + Zh)R a » - (h + 2k 2 )8»R} 



3 v 

2^G a », (22) 



provided that 



fci + 2k 2 + 3k 3 = 6, (23) 

h + 2k 2 = 3. (24) 

For the ADM mass at the spatial infinity, let's use the Schwarzschild metric in Carte- 
sian coordinate 

ds 2 = - 1 1 —J dt 2 + ( 1 + — ^— J (dx 2 + dy 2 + dz 2 ). (25) 

For the energy-momentum four vector 

2^P a = ~Uj"% u , (26) 
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the associated ADM mass energy term is 



M = -i-[/ tH e ^ = 1 ( kl + h ) M , (27) 

which generated one more constraint 

h + k 3 = 2. (28) 
Consider (1231 . fT24"|) and fl28|) . they provide an unique solution 

h = k 2 = k 3 = 1. (29) 

Therefore only the Einstein pseudotensor or other similar ones, such as the Landau- 
Lifshitz pseudotensor, have this property. But when one examines the gravitational 
energy for Landau-Lifshitz or Einstein, they both do not have the good results in 
vacuum, namely only the Bel- Robinson tensor [3j H] . It seems hopeless to have the 
desired result unless one introduces the flat metric tensor rj a p along with g a p. The 
details will be discussed in the next section. 

4 The calculation results 

4.1 The Einstein and Landau-Lifshitz pseudotensors 

The superpotentials for the Einstein eU ^ and Landau-Lifshitz LU a ^ u \ and Bergman- 
Thomson bU 01 ^ classical pseudotensors are 

eU W = -y/^g^r^Z, (30) 
L irW = ^Tg B U^ = gg af3 g^T\J^. (31) 

The pseudotensors can be obtained respectively as follows 

t^ = d v U a ^ v \ t a » = d u U alM - (32) 

The result inside matter at the origin for the Einstein and Bergmann-Thomson pseu- 
dotensors are 

2KtJ{0) = 2Gj{0) = 2kT/(0), (33) 
so their corresponding energy densities are 

£ = -C = -^l = -To = p, (34) 



where p is the mass-energy density. This feature is important because we have to 
fulfill this basic requirement according to the equivalent principle. This means the 
pseudotensor has to match the energy density inside matter at the origin as the metric 
tensor becomes flat. If the zeroth order is faulty it is not appropriate to study the 
second derivatives which give the gravitational energy in vacuum. According to [3], 
the non- vanishing terms of the Einstein and Landau-Lifshitz pseudotensors in vacuum 
are 

eJ = -2r Ara r^ + ^r AffT r^, (35) 
L a/3 = r a Xa (r f3XrT -r Xaf3 )-r Xa a (rP Xa + r aXf3 ) + g ap r XfTT r' 7X ' r . (36) 

The second derivatives of the Einstein pseudotensor [I] in vacuum is 



— gi^Bap^ - S af3fiu ). (37) 



According to [3] , a similar calculation gives the non- vanishing second order contribu- 
tion for the Landau-Lifshitz L a p classical pseudotensor in vacuum as 



,2 1 (- 1 



d^Lafi — - ( 7B a/3lJlu + -Sappy) . (38) 
In order to obtain a positive multiple of B a p^ ui they have made the combination 

This result is good because it guarantees positivity energy for all reference frames. 

Looking at (|37|) . (|38|) and (|39|) . what do they mean physically? One way of getting 
some physical insight is to integrate over a small coordinate sphere at time t = 0, the 
four momentum is 

P )1 = (-E,P t ) = Jt\d 3 x, (40) 

where the energy E should be non-negative. However is not necessarily a proper 
four vector, because it refers to the energy-momentum in a unit sphere of r = constant 
at t = constant and needs not refer to the energy-momentum in a unit sphere of 
r' = constant at t' = constant. These two are not related by a simple Lorentz 
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transformation. Consider the four momentum to second order for the Einstein pseu- 
dotensor in vacuum within a small sphere, as it appears to the laboratory frame 
observer; it is 

1 f 1 



eP, = Yk J ^(45% -SP^zViPx 



540G 

„5 



(7E ab E ab - 3H ab H a \ Ae c ab E ad H\), (41) 



since 



540G 

1 ,„•„• f o ,Q 47T 



I xVd 3 z = | r 2 d 3 x = -^r 5 ^'. (42) 
Similarly for the Landau-Lifshitz pseudotensor 



lPu = — I — ( 7B° uij + -S ^ ) xV(2 3 £ 
M 2k J 18 V m 2 ™ 

~5 



1080G 

^(^^ + 6tf a6 tf afc , 7e c ab E ad H\). (43) 
Combining these two pseudotensors, the properly normalized superpotential becomes 

EL [/a[H = 2 |I^ £[/a [H + B [/a[Hj . (44) 

The gravitational energy-momentum within a small sphere is 

2 n 



r 5 i 



180G 

-XE ab E ab + # ab iT b , 2e c ab i? a(i i/ d 5 ). (45) 



180G 



This result gives the coefficient of the energy- momentum density should be. In other 
words, the coefficient of the Bel-Robinson tensor in vacuum in a holonomic frame is 

fytj = \B a \ v . (46) 
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4.2 The modification of the Einstein pseudotensor 

The Einstein pseudotensor offers a good foundation for the gravitational energy ex- 
pression. One may wonder that if this pseudotensor can be modified so that it has a 
nice positive gravitational energy result. Based on the simple and natural boundary 
conditions, analogy with the Dirichlet or Neumann boundary conditions, one of the 
Chen-Nester four quasilocal expressions can satisfy this requirement [2J. 

According to [2J, there is a two parameters set for the modified Chen-Nester 
quasilocal boundary expressions. There are infinite number of solutions to achieve a 
positive Bel-Robinson tensor. The quasilocal boundary expressions, in compact form, 
are 

2k8 CuC2 {N) = 2nB p (N) + cii N Ar a pAAr] a l3 -C2Ar a pAi N Ari a > 3 

where ci, c 2 G 9ft and h a p := g a p — r) a p. When (ci, c 2 ) = (0, 0), (0, 1), (1, 0) and (1,1), 
this recovers the original Chen-Nester four expressions. From (jUJ), the superpotential 
can be extracted as 

t/ Q [H = E U a [fMu] + c^gh^T^JZ + c 2 ^-gh^Y T ^5 x T Z. (48) 

This superpotential looks like a modification of the Freud superpotential. Note that 
the hT terms do not affect the results inside matter and at spatial infinity, but they 
do affect second order vacuum value. As mentioned before, the pseudotensor can be 
obtained as 

t a » = d v U a ^. (49) 
For all ci and c 2 , inside matter at the origin the result is as desired: 

2KtJ(0) = 2G/(0) = 2kTJ(0), (50) 

and the corresponding energy density is 

£ = ~t ° = -— = -T ° = p. (51) 



The second derivatives in vacuum are 
1 
9 



3U/(0) = - {(4 + Cl - 5c 2 )5 a V - (1 - 2 Cl + c 2 )Sj^ + (ci - 3c 2 )(Y a V + 9T a ^ u ) 

(52) 
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and the gravitational energy-momentum is 



(4 + ci - 5c 2 )5%- - (1 - 2 Cl + c 2 )S° A 



p _ _ _ ) V* °1 ~~ °°2^ M ij - ^ - ^1 T c 2 ;^ /iij I <_j j3_ 

2^7 18 1 +( Cl -3c 2 )(y%+9T%) / 



™ 5 



540G 



2 + ci - 4c 2 )^ / + 5(1 - ci - 2c 2 )7) t0 /} . (53) 



In order to obtain the positive gravitational energy, the coefficients of B^ 1 should be 
positive and that of T^i 1 should vanish. Namely 



or equivalently 



Then rewriting f[5"3"j) 



2 + ci - 4c 2 > 0, 1 - ci - 2c 2 = 0, (54) 



ci > 0, (ci,c 2 ) = (ci, 1 Cl ) . (55) 



P _^1r * 
M ~ 180G ^ 

5 

" T^G {EabEab + HabHab, 2e - abE -d Hd b)- ( 56 ) 

There are an infinite number of solutions since there is a parameter c\ which one can 
tune. Different solutions represent different boundary conditions which is explained 
in [2]. However there is one solution with positive energy and a simple boundary 
condition which is when (ci,c 2 ) = (1,0) in ( 1471) . It reduces to one of the original 
Chen-Nester expressions that they called B C (N). In detail 

2kB c {N) = i N T a p A A-qJ + AT a p A i N r]J, (57) 

the corresponding superpotential is 

= E u a w + ^h x *r\JZ- (58) 

This modification is not just achieving the Bel-Robinson tensor by using an artificial 
combination, it comes from a simple Lagrangian derivation and a Hamiltonian anal- 
ysis beginning from the least action principle. Inside matter at the origin and the 
second derivatives of the small sphere in vacuum results respectively are 

2/^/(0) = 2(3/(0) = 2kTJ(0), (59) 



1 



9u V t a p(0) — ■x(^B a j3 flv + Sapnu + Yot^v + 9T a/ g M „). (60) 

y 
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The result inside matter is good. But the second derivatives of the pseudotensor does 
not give only the Bel-Robinson tensor. However after integration, the gravitational 
energy- momentum turns out to be what we expect: 

1 f 1 



180G 



HOI 



r 



(E ab E ab + H ah H a \ 2e c ab E ad H\), (61) 



180G 



which is the same as f l4"5]) . This result shows that the vector P M is future pointing and 
non-spacelike. 



4.3 The modification of the Landau-Lifshitz pseudotensor 

The superpotentials for the Papapetrou classical pseudotensor is 

pH Wla0\ = -^g™ v nb 5X n , (62) 

or equivalently 

pU *[H = dp [ P H l ^ ][af3] ) 

= BU a ^-^(g^h^V a x J^ + g a ^V T x J^), (63) 

At this point, this superpotential looks like a modification of the Bergmann-Thomson 
or Landau-Lifshitz pseudotensor. Once again, the hT terms do not affect the results 
inside matter and at spatial infinity, they would contribute however to the second 
order vacuum value. As mentioned before, the pseudotensor can be treated as 

t a » = d„U a[H - (64) 
Inside matter at the origin and the second derivatives in vacuum [5] are respectively 

2^/(0) = 2G Q /3 (0) = 2kT/(0), (65) 
9%PaM = l^B afiV -S ai , v -Y a ^-9T af3fiV ). (66) 

The result inside matter is good but the second derivative term does not give the de- 
sired result. However after integration, the gravitational energy-momentum becomes 
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good: 

1 r 1 



p = / _ (ab° S° Y° 9T° ■ r i r j d 3 r 

2k J 9 ^ m m m ' 



180G 



5 



r 



5 



£ afe £ afc + # ab iT 6 , 2e c ab E ad H\). (67) 



180G 



Once again, it is the same as ( )45i) and, as said previously, the vector P M is future 
pointing and non- spacelike. 



5 Conclusion 

The Einstein and Landau-Lifshitz pseudotensors do not give a positive gravitational 
energy in vacuum in terms to second order. In order to obtain the expected result 
of the Bel-Robinson tensor in vacuum, one may consider the modification of these 
two pseudotensors. One of the possible modifications of the Einstein pseudotensor 
turns out to be one of the Chen-Nester quasilocal expressions. For the modification of 
the Landau-Lifshitz pseudotensor, one of the choices is the Papapretou pseudotensor. 
There are two points that need to be emphasized, they both have the same magnitude 
of the B^qi term and they both correspond to natural boundary conditions. 
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